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Abstract. We present a semiclassical approach to the SU (N) Yang-Mills theory whose partition function at nonzero 
temperatures is approximated by a saddle point - an ensemble of an infinite number of interacting dyons of N kinds. The 
ensemble is governed by an exactly solvable 3d quantum field theory, allowing calculation of correlations functions relevant 
to confinement. We show that known criteria of confinement are satisfied in this semiclassical approximation: (i) the average 
Polyakov line is zero below some critical temperature, and nonzero above it, (ii) a quark-antiquark pair has linear rising 
potential energy, (iii) the average spatial Wilson loop falls off exponentially with the area, (iv) N 2 gluons are canceled out 
from the spectrum, (v) the critical deconfinement temperature is in good agreement with lattice data. 

Using the same approximation, we find confinement for the exceptional gauge group G(2) and a first-order deconfinement 
transition, also in agreement with lattice findings. 
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INTRODUCTION 

Quark confinement is one of the most puzzling phenom- 
ena in modern physics. It is widely believed (and sup- 
ported by numerical simulations on a lattice) that if one 
attempts to separate a probe quark from an antiquark, or 
a quark from two other quarks in a proton, a force of 
about 14 tons (called the string tension) pulls it back, 
and this force does not decrease with the separation! Al- 
though long-range forces can be found in condensed mat- 
ter physics (e.g. in the Ising model), it is for the first 
time in history we encounter such a phenomenon in a 
local microscopic theory, and it cries for an explanation. 
However, after 35 years QJJ] of Quantum Chromodynam- 
ics (QCD) not only an ab initio mathematical derivation 
of the "14 tons" is absent but there is still no consensus 
on what is the qualitative mechanism of quark confine- 
ment. 

Today we know for sure that nuclear and subnuclear 
physics is governed by QCD - a local renormalizable 
non-Abelian gauge theory of "colored" quarks interact- 
ing with gluons being analogs of photons. Therefore, the 
problem of confinement is, at least, formulated mathe- 
matically: it must be a property of the quantum Yang- 
Mills (YM) theory. Non-Abelian gauge invariance re- 
quires that gluons interact, even in the absence of quarks. 
In this paper we consider the pure YM theory with no dy- 
namical quarks, also called the "pure glue" theory. How- 
ever, quark sources or probes can be inserted to test cer- 
tain YM correlation functions of interest, in particular we 
shall be interested in the correlation function of Polyakov 
lines and the average of large Wilson loops: these are 
quantities that measure the force between probe quarks. 

Gluons are massless due to gauge invariance, and the 



only freedom a theorist has in a pure glue theory is the 
choice of the gauge group. It is SU (3) for the real world, 
but it is helpful to consider Yang-Mills theories based on 
other Lie groups, in particular SU (N) at any N. (At the 
end of the paper we consider also the exceptional G(2) 
group for reasons explained there.) The scaling of most 
physical observables with N at large N can be found from 
simple Accounting rules, therefore considering arbitrary 
SU(N) gauge groups provides a powerful check. 

The gauge coupling constant a s , the analog of the 
QED coupling a ~ 1/137, is in fact not a constant and 
hence not a parameter that one can choose at will: it 
"runs" as function of the characteristic momentum at 
hand. Owing to asymptotic freedom, it is small at large 
and large at small momenta. As in any other asymptot- 
ically free theory with no explicit scale parameter, the 
"dimensional transmutation" occurs in the YM theory: 
an exponentially large correlation length £ appears, be- 
ing the renormalization-invariant combination of the ul- 
traviolet cutoff (e.g. the lattice spacing a) and the bare 
coupling constant a s (a) given at that cutoff, 

, 1 / 3 1 \ ... a,(a)N 

where A (a) is the so-called 't Hooft bare coupling con- 
stant; it does not depend on N at large N. A has the di- 
mension of mass; it is called the YM scale parameter and 
it actually defines all dimensional quantities in the the- 
ory. The deconfinement phase transition temperature T c 
is proportional to A and the string tension is proportional 
to A 2 by dimensions. All dimensionless quantities are, 
generally speaking, of the order of unity, hence it is a 
strong-coupling problem from the start. It makes the the- 
orist's life hard. 



Unfortunately, QCD will hardly be ever proved to 
be an exactly solvable quantum field theory, even in 
the large N limit. Therefore, one can either do exact 
calculations in a theory that has more symmetries but is 
not our world (e.g. considering supersymmetric versions 
of QCD), or work with QCD but make approximations. 
The first is useful as a theoretical laboratory, the second 
is necessary to understand semi-quantitatively the key 
phenomena, to explain experimental data, and to make 
predictions. 

An approximation is considered to be legitimate if 
there is a systematic way of improving its accuracy. The 
semiclassical approach which we develop below, belongs 
to this category. One chooses a saddle-point classical 
field and then has to take into account quantum fluctu- 
ations about it. Part of the fluctuations are ultra-violet 
and are thus the same as in empty space. Their role is to 
renormalize the bare coupling constant; at this point the 
YM scale parameter (Q~|i emerges. What is left, is a se- 
ries in 't Hooft's running coupling X coming from loop 
expansion in the background of classical configurations. 
The running coupling is evaluated at the maximal mo- 
mentum in the problem, be it temperature or the average 
density of classical configurations. Speaking generally, 
such expansion parameter is of the order of unity, how- 
ever numerically it turns out to be small: A is between 
\ at zero temperature and i near T c . Therefore, in the 
whole range of temperatures within the confining phase 
the semiclassical approximation is expected to yield an 
accuracy of 15-25%, already in the 1-loop approxima- 
tion (provided the saddle point is chosen correctly!) with 
a potential for rapid improvement when higher loops are 
taken into account. We shall see, however, that the actual 
accuracy can be much better than this estimate. It is not a 
too big price to pay if confinement, the most challenging 
riddle in 35 years, is explained in simple terms. 

We consider the pure Yang-Mills theory based on the 
SU(N) gauge group in a broad range of temperatures 
between and T c , the deconfinement phase transition 
temperature. Although the formalism we use is designed 
for nonzero T, we shall see that the physical observables 
we find (such as the string tension) have a finite limit 
when T — ► 0. 

Confinement, as we understand it today and learn from 
lattice experiments with a pure glue theory, has in fact 
many facets, and all have to be explained. For exam- 
ple, in a general SU(N) group one can consider "quarks" 
in various irreducible representations. From the confine- 
ment viewpoint all representations are characterized by 
the phase it acquires under the gauge transformation 
from the group center. The representation is said to have 
'W-ality = k" if the phase is 2M. Let us formulate math- 
ematically the main confinement requirements that need 



to be satisfied: 

■ the average Polyakov line in any A^-ality nonzero 
representation of the SU(N) group is zero below T c 
and nonzero above it 

• the potential energy of two static color sources 
(defined through the correlation function of two 
Polyakov lines) asymptotically rises linearly with 
the separation; the slope called the string tension de- 
pends only on the A^-ality of the sources 

• the average of the spatial Wilson loop decays expo- 
nentially with the area spanning the contour; at van- 
ishing temperatures the spatial ("magnetic") string 
tension has to coincide with the "electric" one, for 
all representations 

• the mass gap: no massless gluons should be left in 
the spectrum. 

Remarkably, all these requirements are satisfied al- 
ready in a semiclassical approximation if one uses an 
ensemble of dyons as a saddle point in the Yang-Mills 
partition function J2|] . 



YANG-MILLS THEORY AT NONZERO 
TEMPERATURES 

Following Feynman, the Yang-Mills (YM) partition 
function can be written as a functional integral over 
the YM 4-potentials A^(t,x) that are traceless hermitian 
NxN matrices, satisfying periodic boundary conditions 
in imaginary (or Euclidean) time: 

2?=jDA ll (t,x) exp L^JJdtjd 3 xTrF^X(2) 



An\t+—,xj =A )1 (t,x), T = temperature, 

where F^ v — d^A v - <9 V A M — /[A M A V ] is the YM field 
strength. 

An important variable is the Polyakov loop; it is the 
path-ordered exponent in time direction, hence it can 
depend only on the space point x: 

L(x) = ^exp (^i Jj dtA 4 (t,x)^j , LeSU{N). 

(3) 

Under space-dependent gauge transformations it trans- 
forms as L — > U~ l LU. The eigenvalues of L(x) are gauge 
invariant; we parameterize them as 



L = diag(e 2ni ^,e 27Ci ^,...,e 



(4) 



Hi + . . . + jUjv = 0, and assume that the phases of these 
eigenvalues are ordered: < < . . . < jj. N < = 



FIGURE 1. The perturbative potential energy as function of the Polyakov line for the SU (2) (left) and SU (3) (right) groups. It 
has minima where the Polyakov loop is one of the N elements of the center Zn and is maximal at the "confining" holonomy. 



+ We shall call the set of N phases {/!,„} the "holon- 
omy" for short. Apparently, shifting jx's by integers does 
not change the eigenvalues, hence all quantities have to 
be periodic in all /x's with a period equal to unity. 

The holonomy is said to be "trivial" if L belongs to one 
of the N elements of the group center Zn- For example, 
in SU (3) the three trivial holonomies are 



Mi = M2 = M3 = 
2 1 



L = U 



Mi 
Mi 



-3,M2=3,M3-3 



1 



1 



-3,M2 = -3-,M3 = 3 



L = eT-l 3 , 



L = e~^~h 



Trivial holonomy corresponds to equal /x's, modulo 
unity. Out of all possible combinations of fi's a distin- 
guished role is played by equidistant /i's corresponding 
toTr/. 0: 



M, 



conf . 



11m 

~2~2N + N' 



(5) 



For example, in 51/(3) it is 



1 1 

Mi =~3>M2 = 0,^3=3- 



L = diag (e-^,l,eT-J 



(6) 

We shall call it "most non-trivial" or "confining" holon- 
omy as it corresponds to TrL = which is the 1 st con- 
finement requirement. 

Immediately, an interesting question arises: Imagine 
we take the YM partition function (O and integrate 
out all degrees of freedom except the eigenvalues {/!„,} 
of the Polyakov loop L(x) which, in addition, we take 
slowly varying in space. What is the effective action foi- 
l's? What set of fi's is preferred dynamically by the YM 
system of fields? 

First of all, one can address this question in perturba- 
tion theory: the result for the potential energy as function 



of/x's is SI 



ppert . 



(27T 



,2 r 3 



N 

£(m» 



-Mn) 2 [l-(Mm-Mn)] : 



(7) 



mod 1 



It is proportional to T 3 (by dimensions) and has exactly 
N zero minima when all /x's are equal modulo unity, see 
Fig. 1 . Hence, f pert says that at least at high temperatures 
the system prefers one of the N trivial holonomies corre- 
sponding to the Polyakov loop being one of the N ele- 
ments of the center Z^. However, terms with gradients of 
jU's in the effective action become negative near "trivial" 
holonomy, signalling its instability even in perturbation 
theory H]. 

It is interesting that in the supersymmetric jV = 1 ver- 
sion of the YM theory (where in addition to gluons there 
are spin-^ gluinos in the adjoint representation) the per- 
turbative potential energy Q is absent in all orders ow- 
ing to fermion-boson cancelation, but the nonperturba- 
tive potential energy is nonzero. Moreover, it is known 
exactly as function of /z's Ja]: it has a single minimum 
at precisely the "most non-trivial" or "confining" holon- 
omy (O. The result can be traced to the semiclassical 
contribution of dyons, which turns out to be exact owing 
to supersymmetry. 

In the non-supersymmetric pure YM theory, the dyon- 
induced contribution cannot be computed exactly but 
only in the semiclassical approximation (this is what the 
paper is about), and the perturbative contribution (O is 
present, too. We shall show below that a semiclassical 
configuration - an ensemble of dyons with quantum fluc- 
tuations about it - generates a nonperturbative free en- 
ergy shown in Fig. 2. It has the opposite behavior of the 
perturbative one, having the minimum at the equidistant 
(confining) values of the fi's. There is a fight between the 
perturbative and nonperturbative contributions to the free 
energy [7]. Since the perturbative contribution to the free 
energy is ~ T 4 with respect to the nonperturbative one, it 




FIGURE 2. The dyon-induced nonperturbative potential energy as function of the Polyakov line for the SU (2) (left) and SU (3) 
(right) groups. Contrary to the perturbative potential energy, it has a single and non-degenerate minimum at the confining holonomy 
corresponding to TrL = 0. 



certainly wins when temperatures are high enough, and 
the system is then forced into one of the N vacua thus 
breaking spontaneously the symmetry. At low tem- 
peratures the nonperturbative contribution prevails forc- 
ing the system into the confining vacuum. At a critical T c 
there is a confinement-deconfinement phase transition. It 
turns out to be of the second order for N=2 but first order 
for N = 3 and higher, in agreement with lattice findings. 



DYON SADDLE POINTS 

Dyons or Bogomolny-Prasad-Sommerfield (BPS) 
monopoles [8] are (anti) self-dual solutions of the 
nonlinear Maxwell equations, D"^F^ V = 0. In SU(N) 
there are exactly N kinds of 'fundamental' dyons with 
Coulomb asymptotics for both electric and magnetic 
fields (hence the term "dyon"): 



±E = B 



2|x 



diag(l,-l,0,. 
diag(0,l,-l,. 

diag(0,0,0,.... 
diag(- 1,0,0,. 



.,0,0) 
.,0,0) 

1-1) 

.,0,1) 



(8) 



Dyon solutions are labeled by the holonomy or the set of 
jlTs at spatial infinity: 



A 4 (|x| 



2jrrdiag(/ii,/i 2) . • 



(9) 



The explicit expressions for the solutions in various 
gauges can be found e.g. in the Appendix of Ref. |9[]. In- 
side the cores which are of the size ~ l/(Tv m ), the fields 
are large, nonlinearity is essential. The action density is 
time-independent everywhere and is proportional to the 
temperature. Isolated dyons are thus 3d objects but with 
finite action independent of temperature: 

27T r-> 

Sdyon = V,„, V m = Hm+1 — jU m , £^V m = l 1 (10) 

a s m 



(here /i/v+i = Hi + 1). The full action of all N kinds 
of well-separated dyons together is that of one standard 
instanton: 5i nst = 2n/a s . 

In the semiclassical approach, one has first of all to 
find the statistical weight with which a given classical 
configuration enters the partition function. It is given 
by exp(— Action), times the determinant -1 / 2 from small 
quantum oscillations about the saddle point. For an iso- 
lated dyon as a saddle-point configuration, this factor di- 
verges linearly in the infrared region owing to the slow 
Coulomb decrease of the dyon field ([8j. It means that iso- 
lated dyons are not acceptable as saddle points: they have 
zero weight, despite finite classical action. However, one 
may look for classical solutions that are superpositions 
of N fundamental dyons, with zero net magnetic charge. 
The small-oscillation determinant must be infrared-finite 
for such classical solutions, if they exist. 



INSTANTONS WITH NON-TRIVIAL 
HOLONOMY 

The needed classical solution has been found a decade 
ago by Kraan and van Baal IToll and independently by 
Lee and Lu 111 111 , see also II 1 211 . We shall call them for 
short the "KvBLL instantons"; an alternative name is 
"calorons with nontrivial holonomy". The solution was 
first found for the SU (2) group but soon generalized to 
an arbitrary SU(N) 01311 . see d 1 411 for a review. 

The general solution A^ vBLL depends on Euclidean 
time t and space x and is parameterized by 3N positions 
of N kinds of 'constituent' dyons in space xj, . . . ,Xjy and 

their t/(l) phases , , All in all, there are AN 

collective coordinates characterizing the solution (called 
the moduli space), of which the action 5; nst = 2n/a s 
is in fact independent, as it should be for a general 
solution with a unity topological charge. The solution 



FIGURE 3. Action density inside the SU (3) KvBLL instanton as function of time and one space coordinate, for large (left), 
intermediate (middle) and small (right) separations between the three constituent dyons. 



also depends explicitly on temperature T and on the 
holonomy /Ji,.. . ,Hn- 

Ajf BLL =A a ^(t,x;xi,...,x N ,Yi,...,Y N ; T,n h ...,y. N ). 

(11) 

The solution is a relatively simple expression given by 
elementary functions. If the holonomy is trivial (all /x's 
are equal modulo unity) the expression takes the form of 
the strictly periodic 0(3) symmetric caloron [15] reduc- 
ing further to the standard 0(4) symmetric BPST instan- 
ton fl6ll in the T — ► limit. At small temperatures but 
arbitrary holonomy, the KvBLL instanton also has only 
a small G{T) difference with the standard instanton. 

One can plot the action density of the KvBLL instan- 
ton in various corners of the parameter (moduli) space, 
see Fig. 3. 

When all dyons are far apart one observes N static 
(i.e. time-independent) objects, the isolated dyons. As 
they merge, the configuration is not static anymore, it 
becomes a process in time. In the limiting case of a 
complete merger, the configuration becomes a Ad lump 
resembling the standard instanton. The full (integrated) 
action is exactly the same 5; nst = 2n/a s for any choice 
of the dyon separations. It means that classically dyons 
do not interact. However, they do experience a peculiar 
interaction at the quantum level to which we proceed. 

QUANTUM WEIGHT OF MANY DYONS 

Remarkably, the small-oscillation determinant about a 
single KvBLL SU(N) instanton made of N different- 
kind dyons can be computed exactly lfl7l[l8ll . With this 
experience, the quantum weight of an arbitrary number 
of dyons of N kinds has been suggested in Ref. |2|. In 
the YM partition function, there are saddle points cor- 
responding to any set of K m dyons. In the thermody- 
namic limit V — > °° one needs to take a saddle point with 
ffiV) dyons. Let K m be the number of dyons of kind 
m (m = 1 . . . AO and let x„„ be the coordinate of the / th 



dyon of kind m (i = 1 . . .K m ). In the semiclassical ap- 
proximation the YM partition function is approximated 
by the partition function of a grand canonical ensemble 
of Ki + K 2 + . . -+K N dyons, 

S = E v i 1 v i ft fl / Vdetg(x mi ), 

K i ...K N K V-- K ^- m =\ i=\ J 

(12) 

where g(x mi ) is a A(K\ + . . . + K N ) x A(Ki + . . . + K N ) 
metric tensor of the dyons' moduli space, composed by 
the overlaps of zero modes of individual dyons, and / is 
the fugacity, 

N 2 A 4 

The bare 't Hooft coupling constant X is renormalized 
and starts to "run" only at the 2-loop level not considered 
here. Eventually, its argument will be the largest scale 
in the problem, be it the temperature or the equilibrium 
density of dyons. 

It is not difficult to find the metric tensor g(x m i) for 
well-separated dyons. In this case the four zero modes 

0fj = 1,2,3,4) of individual dyons are given by the 

Ik) 

components of the field strength: 0^ ' = F^ K . The zero 
modes for the m th kind of dyon are normalized to its ac- 
tion, /Tr^ K) 4 A) - 8 Kl V,„ (see Eq. OB) and hence de- 
pend on the holonomy. Since the field strengths decay 
as 1/r 2 (see Eq. ([8j) the overlaps between zero modes 
are Coulomb-like, and only those that are nearest neigh- 
bors in m do interact. In fact, the diagonal components of 
the metric tensor also acquire Coulomb-like corrections 
since the action of individual dyons is actually normal- 
ized to its asymptotic field A4 that gets Coulomb correc- 
tions from other dyons. 

As a result, we obtain the A(K X + . . . + K N ) x A(K\ + 
. . . +Kn) metric tensor g(x m i) with Coulomb interactions 
as entries, and the v m 's on the diagonal. It turns out that 
its determinant is a square of the determinant of a related 



matrix, \/detg = detG where G is a (K\ + . . . +Kn) x 
(K\ + .. .+K N ) matrix: 



Anv,„ + > 



(14) 



$m,n— 1 




^ I x mi x m j | 



where x,„, is the coordinate of the z th dyon of kind m. The 
matrix G has the following nice properties: 

• symmetry: G mi . nj = G„j. mi 

• overall "neutrality": the sum of Coulomb interac- 
tions in non-diagonal entries cancel those on the di- 
agonal: "L„jG mi .„j = 47TV,„ 

• identity loss: dyons of the same kind are indis- 
tinguishable, meaning mathematically that detG is 
symmetric under permutation of any pair of dyons 
(£«-> jf) of the same kind m. Dyons do not 'know' to 
which instanton they belong to 

• attraction/repulsion: if one decreases the separation 
between same-kind dyons or increases the sepa- 
ration between different-kind dyons, the detG de- 
creases. It means that same-kind dyons repulse 
each other whereas different-kind dyons attract each 
other. The det G measure favors formation of neutral 
clusters with N different kinds of dyons 

• factorization: in the geometry when dyons fall into 
K well separated neutral clusters of N dyons of 
different kinds, detG factorizes into a product of 
exact integration measures for K KvBLL instan- 
tons Oil |2(| valid for any separations between 
different-kind dyons, including their strong overlap 

• last but not least, the metric g corresponding to G is 
hyper-Kahler, as it should be for the moduli space of 
a self-dual classical field |2l!]. In fact, it is a severe 
restriction on the metric. 

An overall constant factor depending on the holon- 
omy and temperature, exp(— P pert V), is understood in 
Eq. (fT2l . where f pert is the perturbative gluon loop (|7]i 
in the background of a constant field A4 (0. This factor 
arises from the non-zero modes in the fluctuation deter- 
minant about dyons and is necessarily present as most of 
the 3d space outside the dyons' core is just a constant A4 
background. Indeed the calculation lfl7[fl8ll exhibits this 
factor which is the only one proportional to the 3-volume 
V. 

The ensemble defined by a determinant of a matrix 
whose dimension is the number of particles, is not a 
usual one. More customary, the interaction is given by 



the Boltzmann factor exp ( — f/i nt (xi , ...)). Of course, one 
can always present the determinant in that way using 
the identity detG = exp(TrlogG) = exp(— Ui nt ) but the 
interactions will then include three-, four-, five-... body 
forces. At the same time, it is precisely the determinant 
form of the interaction that makes the statistical physics 
of dyons an exactly solvable problem. 



STATISTICAL PHYSICS OF DYONS AS 
A QUANTUM FIELD THEORY 

It is possible to present the grand canonical ensemble of 
dyons, governed by the interaction ( TBI as an equivalent 
3d quantum field theory. This will enable us to compute 
various correlation functions of interest. 

To proceed to the quantum field theory description we 
use two mathematical tricks. 



1. "Fermionization" (Berezin [22]). It is helpful to ex- 
ponentiate the Coulomb interactions rather than keeping 
them in detG. To that end one presents the determinant 
of a matrix as an integral over a finite number of anti- 
commuting Grassmann variables ({ \\f\ = 8abY 

det(G AB )= jY\d\i/ld\i/ A exp(y/lG A By'B) ■ 



Now we have the two-body Coulomb interactions in 
the exponent and it is possible to use the second trick. 

2. "Bosonization" (Polyakov [23]). One can present 
the Coulomb interactions in the exponent with the help 
of a Gaussian integral over an auxiliary field (j): 



ex P f E u~~ 1 1 = / D § ex P 



-djtpdjij) + ptj) 



dx 



16n 

P =Egm<5(x-X m ) 



4tz 



After applying the first trick the "charges" Q m become 
Grassmann variables but after applying the second one, 
it becomes easy to integrate them out since the square 
of a Grassmann variable is zero. In fact one needs 2N 
boson fields v,„, w„, to reproduce diagonal elements of G 
and 2N anticommuting ("ghost") fields XnnXm to present 
the non-diagonal elements. The chain of identities is 
accomplished in Ref. [2] and the result for the partition 
function for the dyon ensemble ( TT2l is, identically, a path 
integral defining a quantum field theory in 3 dimensions: 

2? = J Dx^DxDvDw exp Jd 3 x^ {diX^iXm 



diV m diW m )+f 
d 2 ^ 



(-4^ m + v m ) 



<9w„, 



<9w m <9w„ 



w m -w m+ i 



(15) 



The fields v,„ have the meaning of the asymptotic 
Abelian electric potentials of dyons, 



A m4 (x)/T = 27Tjit m -Av m (x), E m = VA„ 



while w m have the meaning of the dual (or magnetic) 
Abelian potentials. Note that the kinetic energy for the 
v m ,w m fields has only the mixing term <9,v m <9,w„, which 
is nothing but the Abelian duality transformation E • B. 
The function J^"(w) in (T5[ where one assumes a cyclic 
summation over ;«, is known as the periodic (or affine) 
Toda lattice. 

Although the Lagrangian in Eq. Sl5[ describes a 
highly nonlinear interacting quantum field theory, it is 
in fact exactly solvable! To prove it, one observes that 
the fields v m enter the Lagrangian only linearly, there- 
fore one can integrate them out. It leads to a functional 
5-function: 



Dv„ 



An 



(17) 



This 5-function restricts possible fields w m over which 
one still has to integrate in eq. ( TT3T >. Let w„, be a solution 
to the argument of the 5-function. Integrating over small 
fluctuations about w gives the Jacobian 



Jac = det 



An 



d 2 8 m „ + f 



<9w m <9w„ 



(18) 



Remarkably, exactly the same functional determinant but 
in the numerator arises from integrating over the ghost 
fields, for any background w. Therefore, all quantum 
corrections cancel exactly between the boson and ghost 
fields (a characteristic feature of supersymmetry), and 
the ensemble of dyons is basically governed by a clas- 
sical field theory. 

To find the ground state we examine the fields' poten- 
tial energy being —4nfpL m d,'^/dvi m which we prefer to 
write restoring v„, = ji nl+ 1 — ji m and & as 



& ■ 



-47T/V£v m e 



w,„-w,„ +1 



(19) 



(the volume factor arises for constant fields w,„). One has 
first to find the stationary point in w m for a given set of 
V m 's. It leads to the equations 







whose solution is 



(V\V2V3...Vn) n 
Vl 



etc. 



(20) 



Putting it back into eq. ( [T9| > we obtain 
@> = -AxfVN{ViV 2 ...V N )T> , vi + ... + Vjv= 1. (21) 



(16) -phe minimum equal ^ m in = — AnfV is achieved at V] = 



. . . = Vjv = jr, that is at equidistant, confining value of 
the holonomy, cf. Eq. ©. We have also proven that the 
result is exact, as all potential quantum corrections cancel 
in the partition function ( [T5| >. 

Given this cancelation, the key finding - that the dyon- 
induced free energy has the minimum at the confining 
value of holonomy - is trivial. If all Coulomb interactions 
cancel after integration over dyons' positions, the weight 
of a many-dyon configuration is the same as if they 
were infinitely dilute (although they are not). Then the 
weight, what concerns the holonomy, is proportional to 
the product of diagonal matrix elements of G in the 
dilute limit, that is to the normalization integrals for dyon 
zero modes, that is to the product of the dyon actions 
~ V m where V m = ^i m+i - \i m and V# = \l\ + 1 - [i N 
such that V1 + V2 + ■■■ + Vn ~ 1 • The sum of all N kinds 
of dyons' actions is fixed and equal to the instanton 
action, however, it is the product of actions that defines 
the weight. The product is maximal when all actions 
are equal, hence the equidistant or confining /x's are 
statistically preferred. Thus, the average Polyakov line 
is zero, < TrL >= 0. 



HEAVY QUARK POTENTIAL 

The field-theoretic representation of the dyon ensemble 
enables one to compute various YM correlation functions 
in the semiclassical approximation. The key observables 
relevant to confinement are the correlation function of 
two Polyakov lines (defining the heavy quark potential), 
and the average of large Wilson loops. A detailed calcu- 
lation of these quantities is performed in Ref. 120; here 
we only present the results and discuss the meaning. 



Af-ality and ^-strings 

From the viewpoint of confinement, all irreducible 
representations of the SU (N) group fall into jV classes: 
those that appear in the direct product of any number 
of adjoint representations, and those that appear in the 
direct product of any number of adjoint representations 
with the irreducible representation being the rank-£ anti- 
symmetric tensor, k = 1 , . . . ,N— 1 . 'W-ality" is said to be 



zero in the first case and equal to k in the second. A^-ality- 
zero representations transform trivially under the center 
of the group Z;v; the rest acquire a phase 2nk/N. 

One expects that there is no asymptotic linear potential 
between static color sources in the adjoint representation 
as such sources are screened by gluons. If a representa- 
tion is found in a direct product of some number of ad- 
joint representations and a rank-fc antisymmetric repre- 
sentation, the adjoint ones "cancel out" as they can be all 
screened by an appropriate number of gluons. Therefore, 
from the confinement viewpoint all A^-ality = k represen- 
tations are equivalent and there are only N—l string ten- 
sions OkN being the coefficients in the asymptotic linear 
potential for sources in the antisymmetric rank-A; repre- 
sentation. They are called "fc-strings". 

The value k — 1 corresponds to the fundamental rep- 
resentation whereas k = N— 1 corresponds to the repre- 
sentation conjugate to the fundamental [quarks and anti- 
quarks]. In general, the mnk-(N—k) antisymmetric rep- 
resentation is conjugate to the rank-fc one; it has the same 
dimension and the same string tension, a k ^ = <^N-k,N- 

The behavior of Cfyjv as function of k and N is 
an important issue as it discriminates between vari- 
ous confinement mechanisms. On general Accounting 
grounds one can only infer that at large N and k <C 
N, <7 tN /a LN = {k/N)(l + ff{\/N 2 )). Important, there 
should be no 0{N~ l ) correction l24ll . A popular version 
called "Casimir scaling", according to which the string 
tension is proportional to the Casimir operator for a given 
representation (it stems from an idea that confinement is 
somehow related to the modification of a one-gluon ex- 
change at large distances), does not satisfy this restric- 
tion. 



Correlation function of Polyakov lines 

To find the potential energy V^.n of static "quark" and 
"antiquark" transforming according to the antisymmetric 
rank-fc representation, one has to consider the correlation 
of Polyakov lines in the appropriate representation: 



In the general antisymmetric rank-£ representation 



TrL^zi) TrLl Af (z 2 )) =const. exp 



Vk,jv(zi-Z2) 



(22) 

Far away from dyons' cores the field is Abelian and 
in the field-theoretic language of Eq. ( fT3T > is given by 
Eq. ( fT6] l. Therefore, the Polyakov line in the fundamental 
representation is 



TrL(z) = £ Z m , Z,„ = exp ( 2nifi, n -^-v, n (z) ) . (23) 
m=l V Z 



N 

TrL^z) = £ 

m\ <m2<...<m/ l 



where cyclic summation from 1 to N is assumed. 

The average d22b can be computed from the quantum 
field theory ( [15) . Inserting the two Polyakov lines (|24|) 
into Eq. (fT~5b we observe that the Abelian electric poten- 
tial v m enters linearly in the exponent as before. There- 
fore, it can be integrated out, leading to a 5 -function for 
the dual field w,„, which is now shifted by the source (cf. 
Eq. (O): 



- 5 (x-zi )($„,„, + . . . + 8, nmk ) 



8(x-2. 2 )(8 m „ 1 +■■■ + 8„ mk ) 



One has to find the dual field w m (x) nullifying the argu- 
ment of this 5-function, plug it into the action 



exp 



(25) 



and sum over all sets {mi < W2 < ... < 
m k}, {n\ < «2 < ■■• < n k] with the weight 
exp(27r/(mi + . . .+ntk — n\ — . . .— nuj/N). The Jaco- 
bian from resolving the 5 -function again cancels exactly 
with the determinant arising from ghosts. Therefore, 
the calculation of the correlator (122) . sketched above, is 
exact. 

At large separations between the sources \z\ — wq\, the 
fields w m resolving the 5-function are small and one can 
expand the Toda chain: 



1 



W„ 



v mn vy n i 



(26) 



where 



Ji ■ 



/ 2 -1 ... 
' -1 2 -1 ... 

0-12-1 



V -1 







-1 ^ 



... 

-i 2 J 



(27) 



As apparent from Eq. d26l ). the eigenvalues of jM deter- 
mine the spectrum of the dual fields w,„. There is one 
zero eigenvalue which decouples from everywhere, and 
N—l nonzero eigenvalues 



= 2 sin — 
N 



l,...,N-l. 



(28) 



Certain orthogonality relation imposes the selection 
rule: the asymptotics of the correlation function of two 
Polyakov lines in the antisymmetric rank-A: representa- 
tion is determined by precisely the k th eigenvalue. We 
obtain H 

(TrLwfa) TrL^za)) (29) 
" 2 = °° const, exp (^-|z! - z 2 |Mv / 3#W) 
where M is the 'dual photon' mass, 



Eq. (S3 

may look contradictory as we first use = 
curlA; and then Z?,- = <9,A4. Actually there is no contra- 
diction as the last equation is true up to Dirac string 
singularities which carry away the magnetic flux. If the 
Dirac string pierces the surface spanning the loop it gives 
a quantized contribution exp (2 %i -integer) = 1; one can 
also use the gauge freedom to direct Dirac strings parallel 
to the loop surface in which case there is no contribution 
from the Dirac strings at all. 

Let us take a fiat Wilson loop lying in the (xy) plane at 
Z=0. Then eq. ( 1331 is continued as 



M 



NA Z 



4g/ = 

T 2kXT 



@(N). 



(30) 



Comparing it with the definition of the heavy quark 
potential d22l) we find that there is an asymptotically 
linear potential between static "quarks" in any A^-ality 
nonzero representation, with the ^-string tension 



Gk,N 



%k A N %k 

2MT sin — = — sin — . 

N X K N 

(31) 

This is the so-called 'sine regime': it has been found 
before in certain supersymmetric theories H25I1 . Lattice 
simulations 1 26] support this regime, whereas another 
lattice study M27I1 gives somewhat smaller values but 
within two standard deviations from the values following 
from eq. (I3TI 1. 

We see that at large N and k <C N, ow/ai w = 
(k/N)( 1 + ^(1 /N 2 )), as it should be on general grounds, 
and that all fc-string tensions have a finite limit at zero 
temperature. 



AREA LAW FOR LARGE WILSON 
LOOPS 

The magnetic field of dyons beyond their cores is 
Abelian and is a superposition of the Abelian fields of 
individual dyons. For large Wilson loops we are inter- 
ested in, it is this superposition field of a large number 
of dyons that contributes most as they have a slowly de- 
creasing 1 /|x— x,-| asymptotics, hence the use of the field 
outside the cores is justified. Owing to self-duality, 

{B, MU - [diA 4 (x)] m „ = - 1 8 mn d,v m (x) , (32) 

cf. eq. ( fT6b . Since A,- is Abelian beyond the cores, one 
can use the Stokes theorem for the spatial Wilson loop: 



W 



Tr exp i fAjdx' = Tr exp i / B, d 2 & 



£exp ( -i- / d 2 a'di\„ 



(33) 



W = £exp ( -i— / d 3 xd z v,„8(z) 

2- Jx.vGArea 



= L ex P 'T / d\v m d z S{z)) . (34) 

It means that the average of the Wilson loop in the dyon 
ensemble is given by the partition function ( fl3T > with the 



E ex P 



d^x v, 



dS(z) 
dz 



9(x,y e Area)j 



where 9(x,y £ Area) is a step function equal to unity if 
x,y belong to the area inside the loop and zero otherwise. 

As in the case of the Polyakov lines the presence of the 
Wilson loop shifts the argument of the 5-function arising 
from the integration over the v m variables, and the ghost 
determinant cancels exactly the Jacobian from the fluc- 
tuations of w m 's, therefore the classical-field calculation 
is exact. 

One has to solve the non-linear Toda equations on 
w m 's with a source along the surface of the loop, 



-d 2 w m +M 2 ( e w »'- w »'+i - e^-i-w™) 
d3 (z) 

= -IniSmmi —5 — 9{x,y S Area), 
dz 



(35) 



for all m\, plug it into the action (Anf /N)^(w), and 
sum over m\. In order to evaluate the average of the 
Wilson loop in a general antisymmetric rank-fc repre- 
sentation, one has to take the source in eq. d35l > as 
—2ni8'(z) {pmmi +■ • ■+S m ,„ k ) and sum over mi< ...< 
fiik from 1 to N, see eq. (l24l . 

Contrary to the case of the Polyakov lines, one cannot, 
generally speaking, linearize eq. d35T > in w m but has to 
solve the non-linear equations as they are. The Toda 
equations ( f35t with a 8'(z) source in the r.h.s. define 
"pinned soliton" solutions w m (z) that are Id functions 
in the direction transverse to the surface spanning the 
Wilson loop but do not depend on the coordinates x,y 
provided they are taken inside the loop. Beyond that 
surface w,„ = 0. Along the perimeter of the loop, w m 
interpolate between the soliton and zero. For large areas, 



the action (1251 1 is therefore proportional to the area of 
the surface spanning the loop, which gives the famous 
area law for the average Wilson loop. The coefficient in 
the area law, the 'magnetic' string tension, is found from 
integrating the action density of the soliton w m (z) in the 
z direction. 

The exact solutions of Eq. d3~5l for any N and any 
representation k have been found in Ref. |2|], and the 
resulting 'magnetic' string tension turns out to be 



A 2 N . Ttk 



(36) 



which coincides with the 'electric' string tension OTb 
found from the correlators of the Polyakov lines, for all 
fc-strings! 

Several comments are in order here. 

• The 'electric' and 'magnetic' string tensions should 
coincide only in the limit T — ► where the Eu- 
clidean 0(4) symmetry is restored. Both calcula- 
tions have been in fact performed in that limit as 
we have ignored the temperature-dependent pertur- 
bative potential (jT). If it is included, the 'electric' 
and 'magnetic' string tensions split. 

• despite that the theory $15[ is 3-dimensional, with 
the temperature entering just as a parameter in the 
Lagrangian, it "knows" about the restoration of Eu- 
clidean 0(4) symmetry at T — > 0. 

• the 'electric' and 'magnetic' string tensions are 
technically obtained in very different ways: the first 
is related to the mass of the elementary excitation 
of the dual fields w m , whereas the latter is related to 
the mass of the dual field soliton. 



weak until T sa T c where it abruptly rises owing to the 
excitation of many glueballs. 

The nonperturbative free energy corresponding to the 
minimum of the dyon-induced potential energy as func- 
tion of the holonomy (fJTJ is 



yon 



v 



N z A 4 

'U? x 2 



(38) 



It is 0(N 2 ) but temperature-independent. We have dou- 
bled the minimum from eq. (fJTJ keeping in mind that 
there are also anti-dyons and assuming that their interac- 
tions with dyons is not as strong as the interactions be- 
tween dyons and anti-dyons separately, as induced by the 
determinant measure (1141 1. therefore treating dyons and 
anti-dyons as two independent "liquids". (By the same 
logic, the string tension OTb has to be multiplied by y/2 
as due to anti-dyons.) 

Dyons force the system to have the "most nontrivial" 
holonomy ©. For that holonomy, the perturbative poten- 
tial energy (Q is at its maximum equal to 



pert, max 
V 



7T 4 
= — T 4 
45 



N N 2 



(39) 



The full free energy is the sum of the three terms above. 

We see that the leading ff(N 2 ) term in the Stefan- 
Boltzmann law is canceled by the potential energy pre- 
cisely at the confining holonomy point and nowhere else! 
In fact it seems to be the only way how ff(N 2 ) massless 
gluons can be canceled out of the free energy, and the 
main question shifts to why does the system prefer the 
"most nontrivial" holonomy. Dyons seem to answer that 
question. 



CANCELATION OF GLUONS IN THE 
CONFINEMENT PHASE 



DECONFINEMENT PHASE 
TRANSITION 



To prove confinement, it is insufficient to demonstrate 
the area law for large Wilson loops and the zero average 
for the Polyakov line: it must be shown that there are 
no massless gluons left in the spectrum. We give an 
argument that this indeed happens in the dyon vacuum. 

A manifestation of massless gluons in perturbation 
theory is the Stefan-Boltzmann law for the free energy: 



V 



logiTs 



SB 



V 



45 



T 4 (N 2 -l). (37) 



It is proportional to the number of gluons N 2 — l and 
has the T 4 behavior characteristic of massless particles. 
In the confinement phase, neither is permissible: If only 
glueballs are left in the spectrum the free energy must 
be ff(N°) and the temperature dependence must be very 



As the temperature rises, the perturbative free energy 
grows as T 4 and eventually it overcomes the negative 
nonperturbative free energy d38l . At this point, the triv- 
ial holonomy for which both the perturbative and non- 
perturbative free energy are zero, becomes favorable. 
Therefore an estimate of the critical deconfinement tem- 
perature comes from equating the sum of Eq. ( 1381 and 
Eq. ( 1391 to zero, which gives 



T 4 = 



45 N 4 



27T 4 N 4 - 1 X 2 



(40) 



As expected, it is stable in N. A more robust quantity, 
both from the theoretical and lattice viewpoints, is the ra- 
tio T c / y/a where a is the string tension in the fundamen- 
tal representation, since in this ratio the poorly known 





SU{3) 


SU(4) 


SU{6) 


SU (8) 




T c j y/a, theory 


0.6430 


0.6150 


0.5967 


0.5906 


0.5830+ °^ 95 + °^ 06 + ... 


Tc/y/a, lattice 


0.6462(30) 


0.6344(81) 


0.6101(51) 


0.5928(107) 


0.5970(38) + ^2) (fit) 



parameters A and X cancel out: 

T c _ ( 45 n 2 N 2 \ * >°q i/45\' 
v/ff" l^ 4 (jV 4 -l)sin 2 ^y ~^ it\4/ ' 

(41) 

In the Table, we compare the values from Eq. ( f4~TT > to 
those measured in lattice simulations of the pure SU(N) 
gauge theories [28]; there is a surprisingly good agree- 
ment. A detailed study of the thermodynamics of the 
phase transition will be published elsewhere. 

AN EXCEPTIONAL GIFT 

It is illuminating to go beyond the SU (N) gauge groups 
and consider the YM theory based on e.g. the exceptional 
group G(2). It has rank 2 and is similar to the SI/ (3) 
group but has only a trivial center, meaning that only a 
unity element commutes with all other group elements 
(in SU(N) there are N such matrices being -roots of 
unity; they form the group center Zjy). Since confinement 
in SU (N) is often associated with the symmetric distribu- 
tion of the Polyakov line eigenvalues with respect to the 
Zn permutations, and the deconfinement transition is as- 
sociated with the spontaneous breaking of Zn symmetry, 
while neither is a feature of G(2), questions have been 
raised whether there is at all confinement in G(2), leave 
alone the confinement-deconfinement phase transition. 

However, recent numerical simulations of G(2) per- 
formed by three groups [29] showed that there is con- 
finement at low T and a first-order deconfinement transi- 
tion. By confinement we mean here the zero average of 
the Polyakov line in the lowest 7-dimensional representa- 
tion. These findings pose difficulties for the confinement 
scenarios based on the center of the group, e.g. on center 
vortices. 

The dyon scenario works, however, exceptionally well 
it this case, too. We have repeated all the steps described 
in this paper but applied to the G(2) group. There are 
three types of fundamental dyons, like in SU (3), but neu- 
tral clusters are formed by four; overall neutrality re- 
quires that one type of dyons have to enter twice more of- 
ten than the two others. The minimum of the free energy 
for the dyon ensemble lies exactly at the holonomy cor- 



responding to the zero Polyakov line. At a critical tem- 
perature there is a first-order transition to a phase with a 
nonzero average Polyakov line [30]. 

We stress that these results are sensitive to the dynam- 
ics as they are not in the least enforced by symmetry. 
They provide strong support to the dyon scenario of con- 
finement. It should be added that instantons with nontriv- 
ial holonomy and dyons have been directly observed in 
lattice simulations Bill . 

SUMMARY 

What happens in the semiclassical approximation based 
on dyons, can be summarized as follows: 

■ The ensemble of dyons favors dynamically the con- 
fining value of the holonomy. This is almost clear, 
given that the weight is proportional to the product 
of individual actions of fundamental dyons, and it 
is maximal when the actions are equal. Such holon- 
omy corresponds to the zero of the Polyakov line 

■ Dyons form a sort of Coulomb plasma (but an ex- 
actly solvable variant of it) with an appearance of 
the Debye mass both for "electric" and "magnetic" 
(dual) photons. The first gives rise to the exponen- 
tial fall-off of the correlation of two Polyakov lines, 
i.e. to the linear heavy-quark potential, the second 
yields the area law for spatial Wilson loops 

■ &{N 2 ) massless gluons cancel out from the free 
energy, and only massive (string?) excitations are 
left. 

The reason why a semiclassical approximation works 
well for strong interactions (where all dimensionless 
quantities are, generally speaking, of the order of unity) 
is not altogether clear. A possible justification has been 
outlined in the Introduction: After UV renormalization 
is performed about the classical saddle points and the 
scale parameter A appears as the result of the dimen- 
sional transmutation, further quantum corrections to the 
saddle point are a series in the running 't Hooft coupling 
X whose argument is typically the largest scale in the 
theory, in this case max(r,n'/ 4 ) where n is the 4a! den- 
sity of dyons. An estimate shows that the running X is 



between 1 /4 at zero temperature and 1/7 or less at crit- 
ical temperature. Therefore, although these numbers are 
"of the order of unity", in practical terms they indicate 
that high order loop corrections are not too large. Let us 
recall that quite an accurate computation of anomalous 
dimensions in critical phenomena from the £-expansion 
by Fisher and Wilson 1 32] is based on truncating the Tay- 
lor expansion in e at the first couple of terms, where e = 1 
or sometimes 2 01 

Unfortunately, approximations made in Ref. [2] and 
reproduced above are not limited to neglecting higher 
loop corrections. We have (i) ignored dyon interac- 
tions induced by the small oscillation determinant over 
nonzero modes (although we did take into account that it 
renormalizes the gauge coupling giving rise to the scale 
parameter A, and that it leads to the perturbative potential 
energy as function of the holonomy), (ii) neglected the 
interactions of dyons of opposite duality, treating them 
as two noninteracting "liquids", (iii) conjectured a sim- 
ple form of the dyon measure which may be incorrect 
when two same-kind dyons come close. Although cer- 
tain justification for these approximations can be put for- 
ward H it is desirable not to use them at all, and that 
may be possible. 
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